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ABSTRACT 

We present a morphological study of nebulae around Wolf-Rayet (WR) stars using archival narrow-band optical and Wide-field 
Infrared Survey Explorer (WISE) infrared images. The comparison among WISE images in different bands and optical images proves 
to be a very efficient procedure to identify the nebular emission from WR nebulae, and to disentangle it from that of the ISM material 
along the line of sight. In particular, WR nebulae are clearly detected in the WISE W4 band at 22 pm. Analysis of available mid-IR 
Spitzer spectra shows that the emission in this band is dominated by thermal emission from dust spatially coincident with the thin 
nebular shell or most likely with the leading edge of the nebula. The WR nebulae in our sample present different morphologies that we 
classified into well defined WR bubbles (bubble ^-type nebulae), clumpy and/or disrupted shells (clumpy/disrupted C-type nebulae), 
and material mixed with the diffuse medium (mixed Al-type nebulae). The variety of morphologies presented by WR nebulae shows 
a loose correlation with the central star spectral type, implying that the nebular and stellar evolutions are not simple and may proceed 
according to different sequences and time-lapses. We report the discovery of an obscured shell around WR 35 only detected in the 
infrared. 
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1. Introduction 


Massive stars (Mi > 30 Mq) end up their lives as Wolf- 
Rayet (WR) stars. Before entering the WR phase, these stars 
evolve through the red or yellow supergiant (RSG or YSG) or 
luminous blue variable (LBV) phases, depending on the initial 
stellar mass, and eject their envelopes via copious slow winds 
expanding at 10-100 km s“^. When a fast stellar wind (1000- 
2,000 km s“^) develops during the WR phase, it sweeps up the 
previously ejected slow RSG/LBV wind material to form a WR 
nebula. The swept-up circumstellar material is photoionized by 
the central star and has temperatures of -10"^ K. 


WR nebulae present different morphological features such 
as bubbles, blowouts, clumps, filaments, and diffuse emis¬ 
sion. These nebulae are observable at all frequencies of the 
electromagnetic spectrum: radi o (e.g. , Arnal & Cappal. 1 19961: 


electromagnetic spectrum: radi o (e.g. , Arnal & Gappal. llWol: 
lArnal et all 119991: l^appaetaLl l200^ 20081 l2009l). infrared 
(IR; e.g. . Ivan Buren & McCravl Il988l: iGvaramadze et al 


ElagevetaLl 120111: IStringM 


201 Oat Mauerhan et al.L l20id: Wachter et al.L l20 10 : 


2012: IWachteretaT 


ow et al 


201 Ih. optical (e.g.. IChuL 119821: iTreffers & ChuL Il982a : 


Chu et al.L 119831: iGmendl et aU I 2 OOOI: IStock & BarlowLl^ : 


Eernandez-Martm et al.L l20l'^ and X rays (e.g.. iBochkarev 
19881: IWriggeetal.1 ll994~iWriggeL IT^: IChu et al.l 12003 


Wrigge & ^ndkei{2002nWrigge et al.Ll2005l:IZhekov & Park . 


Wrigge & WenQkeillzU(Jzi:IWr 

201l]:[Toala et al.Ll2012Ll20ll) 


IChul (Il98ll) started a series of works proposing a coher¬ 
ent morphological classification of nebulae associated with WR 
stars as: 7? - radiatively excite d H n region s, E - stellar ejecta, and 
W - wind-blown bubbles (see IChuL 120031 for an updated review 
of the morphology of WR nebulae). The R-type nebulae present 
optical emission lines with widths comparable to those in ordi¬ 
nary H II regions, and are further divided into amorphous H ii 
regions (Rf) and shell-structured H ii regions (R^). The £’-type 


nebulae, characterized by their highly clumpy appearance and 
irregular velocity field, were proposed to form through sudden 
episodes of mass ejection. Einally, the W-type nebulae present 
a thin sheet of gas and filaments curving around a WR star 
which is close to the geometric center of the nebula or offset 
toward the brightest rim. W-type nebulae are predicted by nu¬ 
merical simul ations of the circumstellar rn ed ium (CSM) around 
WR s t ars (e.g..lGarcia-Segura et al.LlT9963IE : lErever et al.l.l2003L 


20061: iToala & Arthuii l2()lll) . Using this classification scheme. 
Stock & Barlowl (l2010l) extended and revised the morphologies 
of a sample of southern WR nebulae. 


IChul(Il98ll) found a correlation, later confirmed bv IChu et~^ 
(1 19831) . between nebular morphology and spectral type of the 
central WR star: WC stars are mostly associated with R^ nebulae, 
whereas the central stars of W nebulae are mostly early WN stars 
and those of E nebulae are exclusively of spectral type WN8. 
This correlation was interpreted in terms of nebular evolution, 
so that W nebulae , at an early nebular phase, will precede R^ 
nebulae (IChu et al.Lll983l) . 


The classification scheme described above is very demand¬ 
ing observationally. It requires both direct line imaging and spec¬ 
troscopic campaigns to investigate the morphology, kinematics, 
and abundances of the ionized material in order to disentangle 
the emission of the WR nebulae from that of the ISM. In the 
recent years, mid-IR observations have proved very useful to 
study WR nebulae or to discov er new nebulae around evolved 
stars (e.g.. IWachter et al.l l201Ql) . In particular, Spitzer observa¬ 
tions in the 24 pm band have been found to be sensitive to the 
WR nebular emission, whereas the near-IR K and mid-IR 8 pm 
bands trac e mainly the emission of material along the line of 
sight (e.g.. lWachter et al.Ll201Ql:[^ingfellow et al.LT20f^ . Mid- 
IR observations of WR nebulae are thus a promising tool to study 
these objects. 
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Toala et aL: WISE morphological study of WR nebulae 


Here we presen t Wide-field Infrared Survey Explorer (WISE; 
IWright et al.L l201Ql) images of a sample of 31 WR stars and com¬ 
pare them with narrow-band optical images (Section 2). This 
observational approach provides a straightforward method to 
identify nebulae around WR stars, clearly separating its emis¬ 
sion from that of the ISM. The nature of the emission of WR 
nebulae in different WISE bands is investigated using mid- 
IR Spitzer spectra (Section 3). Their morphologies have been 
broadly grouped into three morphological types that can be in¬ 
terpreted in the framework of the evolution of the previously 
ejected dense wind through different evolutive paths (Section 4). 


2. Observations 

We have searched the WISE Preliminary Release Database made 
public in April 2011 for IR im ages of stars in the Vl lth catalogue 
of Galactic Wolf-Rayet stars (Ivan der Huchtl l200ll) and found a 
sample of 31 WR stars with available datfl. To compare their IR 
and optical morphologies, we have selected Ha li ne images of 
these nebulae in the Super COSMOS Sky Survey (iParker et al.L 
l2005b . or STScI Digitized Sky Survey (DSSfrom the Second 
Palomar Observatory Sky Survey (POSS-II). For a few cases, we 
have u sed the [O m] and Ha images presented in iGruendl et al.l 
(|200Q|) that have been kindly provided by Y.-H. Chu and R.A. 
Gruendl. The names and coordinates of the sources in our sam¬ 
ple are listed in Table \T\ together with the distance, height over 
the Galactic Plane, and the spectral type and wind terminal ve¬ 
locity of their central stars. Table [T] also identifies the telescope 
or public database used to acquire the optical images used in this 
paper. Note that for the specific cases of WR7 and WR 136 we 
have also used the existing DSS blue band to complement the 
morphology seen in the DSS red banc0. 

The WISE images of the WR nebulae in Table [T] were down¬ 
loaded from the WISE Image search tool at the NASA/IPAC 
Infrared Science Archive (IRSA). To study the IR morphology of 
these nebulae, we will use the W2 (Ac=4.6pm), W3 (Ac=l2pm), 
and W4 (Ac=22pm) bands. The spatial resolutions for the W2, 
W3, and W4 bands are 6'.'1, 6'.'4, 6'.'5, and 12'.'0, respectively, 
with an astrometric accuracy for bright sources better than O'.' 15. 

2.1. IRS spectra 

To help us interpret the nature of the IR emission of WR neb¬ 
ulae in the WISE W3 and W4 bands, we have searched for 
Spitzer InfraRed Spectrograph (IRS) spectroscopic observations 
in the 10-37 pm range of the WR nebulae in our sample listed 
in Table [T] We found available spectroscopic observations for 
a sample of 15 WR nebulae. Whereas the detailed analysis of 
these spectra will be presented in a subsequent paper (Toala et 
al., in prep.), here we will focus on the nebulae around WR7 
(NGC2359), WR8, and WR31a, because there are available 
high-dispersion spectroscopic observations for these nebulae, as 


^ The mid-IR WISE images of the nebulae around WR 43abc, WR91, 
and WR93 (NGC3603, RCW 122, and NGC6357, respectively) were 
also examined. These stars belong to young stellar clusters, and thus 
their nebulae are influenced not only by the WR star, but by all stars 
in the cluster. Their resulting morphologies is complex and difficult to 
classify. Furthermore, their WISE W3 and W4 band images were satu¬ 
rated at their central regions. Accordingly, we discarded these nebulae 
from our analysis. 

^ https://archive.stsci.edu/cgi-bin/dss_form. 

^ The rest of the WR nebulae do not have their corresponding DSS 
blue band images. 


well as for the background emission from apertures located at 
their periphery. 

The basic calibrated data (BCD) of the IRS observations of 
these WR nebulae were downloaded from the Spitzer Heritage 
Archive. These data are processed with the Spitzer Science 
Center Pipeline Version S 18.18.0, and include high-resolution 
spectroscopic observations (R~600) obtained using the Short- 
High (SH, 9.9-19.6 pm) and Long-High (LH, 18.7-37.2 pm) 
modules. The targets were observed using apertures with sizes 
of 13'.'6x4'.'5 in the SH module and 22'.'3x8'.'9 in the LH mod¬ 
ule. All the IRS data were reduced using the CUbe Builder for 
IRS Spectra Maps {CUBISM). This tool does not only reduce 
the IRS data, but it can also be used to analyze the IRS data 
and to extract one-dimensional spectra. The various processing 
steps followed within CUBISM, including the characterization 
of noise in the data an d the removal of bad pixels, are described 
by ISmith et al.l (l2007l) . We note that the Spitzer IRS spectra of 
WR 8 and WR 31a include the stellar continuum from their WR 
stars, as they were included in the aperture given the small an¬ 
gular size of these nebulae. 

3. Interpreting the IR emission of WR nebuiae 

We present in Figures [T] to [3] and the figures in Appendix 
the optical and WISE IR images of the WR nebulae listed in 
Table [T] The inspection of the colour-composite IR pictures of 
WR nebulae presented in the figures (see Figures [T] to [3]) re¬ 
veals the prevalence of the emission in the WISE W4 band at 
22 pm (red in these colour-composite pictures). The emission 
in this band is mostly coincident with the optical Ha emission 
from ionized nebular material (e.g., WR 16 in Fig.[T]-bottom pan¬ 
els). This is in agreement with Spitzer MIPS 24 pm imaging of 
nebulae around evolved ma s sive s t ars fe.g.. iGyammadze et al.l. 
1200911201 Pal: IWachter et all l2010l: iFlagev et. all l201 Ih . On the 
Other hand, the emission in the W3 band (green in the colour- 
composite pictures), which is sensitive to cold gas or low ion¬ 
ization material, seems to trace mostly ISM gas along the line 
of sight (e.g.. Fig. [3]), with little contribution of emission from 
the WR nebula in some cases (e.g., S 308 around WR6, Fig.[T]). 
Some WR nebulae, however, may show bright emission in this 
band from discrete knots and clumps interior to the nebular shell 
(e.g., NGC6888, Fig. lAl^ . Finally, the W2 band (blue in the 
colour-composite pictures), which is expected to trace the con¬ 
tinuum emissio n from small grains and the background stellar 
component fsee lFlagev et al.L [20111) . shows emission from stars 
in the background. For the most dusty cases, the nebulae around 
WR7, WR18, WR22, and WR23 the W2 band also includes 
some nebular emission. 

The nature of the emission in the WISE W3 and W4 
bands has been further investigated using the Spitzer IRS high- 
dispersion background-subtracted spectra of the nebulae around 
WR7 (NGC2359), WR8, and WR31a presented in Figure [3 
These spectra have been overplotted by the spectral responses of 
the WISE W3 (green) and W4 (red) bands. The spectra of these 
nebulae imply significant contribution to the emission detected 
in the WISE W3 band by a number of emission lines, including 
[S IV] 410.51 pm. He i 410.66 pm. He i 4411.31,12.37,16.21 
pm, [Ar v] 413.10 pm, and [Ne iii] 415.55 pm. In sharp contrast, 
the WISE W4 band is mostly dominated by continuum emission 
probably from thermal dust continuum emission, with weak H i, 
[Fe ii], and [O iv] emission lines. It is interesting to note that the 
spectrum of WR7, extracted from a region known to be over¬ 
abundant in oxygen (iDufouj. [19891) . does not show bright emis¬ 
sion from oxygen lines, but a weak [O iv] 425.59 pm line. 
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Table 1. Stellar and Nebular Parameters of the WR Sample^ 


WR 

Nebula'^ 

R.A. Dec. 

(J2000) 

Galactic Coord. 
l,b 

n 

d 

(kpc) 

z 

(pc) 

spectral 

Type 

Voo 

(kms-b 

Optical 

Observation‘s 

Nebular 

Type 

6 

S308 

06 54 13.05 

-23 55 42.1 

234.76-10.08 

1.50 

-262 

WN4-S 

1700 

CTIO (Her, [0 m]) 

S 

7 

NGC 2359 

07 18 29.13 

-13 13 01.5 

227.75-0.13 

3.67 

-8 

WN4-S 

1600 

SuperCOSMOS/DSS 

S 

8 


07 44 58.22 

-31 54 29.6 

247.07-3.79 

?>A1 

-229 

WN7/WCE+? 

1590 

SuperCOSMOS 

C 

16 

Anon 

09 54 52.91 

-57 43 38.3 

281.08-2.55 

2.37 

-105 

WN8h 

650 

SuperCOSMOS 

S 

18 

NGC 3199 

10 17 02.28 

-57 54 46.9 

283.57-0.97 

2.20 

-37 

WN4-S 

1800 

DSS 

c 

22 

Anon 

10 41 17.52 

-59 40 36.9 

287.17-0.85 

3.24 

-48 

WN7h+09III-V 

1785 

DSS 

M 

23 

Anon 

10 41 38.33 

-58 46 18.8 

286.78-0.03 

3.24 

-2 

WC6 

2340 

DSS 

S 

30 

Anon 

10 51 06.01 

-62 17 01.8 

284.44-2.61 

5.83 

-263 

WC6+06-8 

2100 

SuperCOSMOS 

M 

31a 


10 53 59.66 

-60 26 44.3 

288.94-0.81 

8.0 

-113 

WNllh 

365 

SuperCOSMOS 

S 

35 

Anon 

11 00 22.10 

-61 13 51.0 

289.97-1.19 

17.87 

-371 

WN6h-w 

1100 

SuperCOSMOS 

c 

35b 

Anon 

11 00 02.30 

-60 14 01.0 

289.63-0.24 

2.19 

-9 

WN4 


SuperCOSMOS 

M 

38 

Anon 

11 05 46.52 

-61 13 49.1 

290.57-0.92 

5.83 

-94 

WC4 

3200 

SuperCOSMOS 

M 

40 

RCW58 

11 06 17.21 

-65 30 35.2 

292.31-4.83 

2.26 

-190 

WN8h 

650 

CS (Her, [0 III]) 

C 

52 

Anon 

13 18 28.00 

-58 08 13.6 

306.50+4.54 

1.51 

120 

WC4 

3225 

SuperCOSMOS 

M 

54 

Anon 

13 32 43.79 

-65 01 27.9 

307.27-2.50 

7.53 

-321 

WN5-W 

1500 

SuperCOSMOS 

M 

55 

RCW78 

13 33 30.13 

-62 19 01.2 

307.80+0.16 

6.03 

17 

WN7 (WNE-w) 

1200 

SuperCOSMOS 

C 

68 

G320.5-1.4 

15 18 21.0 

-59 38 10.0 

320.54-1.88 

3.27 

-107 

WC7 

2100 

SuperCOSMOS 

M 

75 

RCW104 

16 24 26.23 

-51 32 06.1 

332.84-1.48 

2.18 

-56 

WN6-S 

2300 

CS (Her, [0 iii]) 

S 

85 

RCW118 

17 14 27.13 

-39 45 47.0 

347.43-0.61 

4.66 

-50 

WN6h-w (WNL) 

1400 

SuperCOSMOS 

c 

86 

RCW 130 

17 18 23.06 

-34 24 30.6 

352.25+1.85 

2.86 

92 

WC7(+B0III-I) 

1855 

DSS 

M 

94 

Anon 

17 33 07.14 

-33 38 23.7 

354.60-0.25 

3.12 

-14 

WN5-W 

1300 

SuperCOSMOS 

M 

95 


17 36 19.76 

-33 26 10.9 

355.13-0.70 

2.09 

-26 

WC9d 

1900 

SuperCOSMOS 

C 

101 

Anon 

17 45 09.10 

-31 50 16.0 

357.47-1.43 

3.18 

-79 

WC8 


SuperCOSMOS 

C 

102 

G2.4+1.4 

17 45 47.00 

-26 10 29.0 

2.38+1.41 

5.56 

137 

W02 

5000 

SuperCOSMOS 

S 

113 

RCW 167 

18 19 07.36 

-11 37 59.2 

18.91+1.75 

1.79 

55 

WC8d+08-9IV 

1700 

DSS 

M 

116 

Anon 

18 27 04.28 

-12 22 52.3 

19.16-0.32 

2.48 

-14 

WN8h 

800 

SuperCOSMOS 

C(M?) 

124 

Ml-67 

19 11 30.88 

+ 16 51 38.2 

50.20+3.31 

3.36 

194 

WN8h 

710 

DSS 

C 

128 

Anon 

19 48 32.20 

+ 18 12 03.7 

55.62-3.79 

9.37 

-619 

WN4(h)-w 

2050 

MLO (Her, [0 iii]) 

s 

131 


20 00 19.12 

+33 15 51.1 

69.90+1.71 

11.78 

352 

WN7h 

1400 

DSS 

c 

134 

Anon 

20 10 14.20 

+36 10 35.1 

73.45+1.55 

1.74 

47 

WN6-S 

1700 

MLO (Her, [0 iii]) 

c 

136 

NGC 6888 

20 12 06.55 

+38 21 17.8 

75.48+2.43 

1.26 

53 

WN6(h)-s 

1600 

DSS (Red, Blue) 

s 


^ All parameters h ave been adopted from Ivan der HuchtI 120011) , but the spectral classification and terminal velocities have been revised according to iHamann ^ inEool) and 
iSanderetaljfMl . 

^The name Anon is given to those nebulae that have not been catalogued as independent objects nebula, as those that are a small portion of a larger nebular region ('see lChulll981l: 
iChu et allfl^ . 

CTIO: Cerro Tololo Inter-American Observatory, SuperCOSMOS: Super COSMOS Sky Survey, DSS: STScI Digitized Sky Survey, CS: Curtis-Schmidt, and MLO: Mount Laguna 
Observatory. 


It is thus concluded that the emission from WR nebulae in 
the WISE W3 band is dominated by H i and forbidden emission 
lines, mostly of [S iv] and [Ne iii], with little contribution to the 
observed emission from dust thermal continuum. Interestingly, 
the spatial correspondence between broad arc-like patches of 
emission in the W3 image of S 308 around WR 6 and diminished 
emission in Ha (Fig. [T]) seems to imply absorption of nebular 
emission by dust on the foreground. As expected, the W3 band 
traces material in the ISM in most cases. Meanwhile, the emis¬ 
sion in the WISE W4 band, tracing mostly the optical nebular 
shells, corresponds to thermal continuum emission, with very 
little contribution from H i, [Fe ii], and [O iv] emission lines. 
These results are not in disagreement with the IRS spectrum of 
the WR bubble MB 3957, which is dominated by emission lines 
of iron, sulfur and neon, but als o includes an addit ional contri¬ 
bution of thermal dust emission (iFlagev et al .11201 ll) . Otherwise, 
the Spitzer MIPS 24 pm images of planetary nebulae include 
important contributions from the [Q iv] 425.8 9 pm and [Ne v] 
424.3 pm emission lines (e.g. IChu et al.[l2009l) . unlike the WISE 
W4 images of WR nebulae. The spatial resolution of the WISE 
W4 band images, ^^1270, does not allow us to determine whether 
the dust in WR nebulae is spatially coincident with the ionized 
material or is found at the leading edge of the optical WR nebu¬ 
lae. 


4. Morphological classification of WR nebulae 

The images of the WR nebulae presented in Figures [H to [3] are 
examples of the variety of morphologies from complete shells, 
arcs (incomplete shells), clumps, and filaments to diffuse, fea¬ 
tureless emission. A close inspection to all WR nebulae in our 
sample and the comparison between optical and IR morpholo¬ 
gies has led us to classify them into three broad morphological 
types: 

- ^ - WR Bubble. 

These nebulae present a thin shell or bubble both in optical 
and IR, mostly in the W4 band. Examples of these nebulae 
are those around WR 6 (S 308) and WR 16 (see Fi^. 

- C - Clumpy/Disrupted WR bubbles. 

These nebulae present clumpy Ha and IR images and/or in¬ 
complete arcs or shells. The W4 images reveal for a signifi¬ 
cant number of cases that the optical and IR clumps are spa¬ 
tially coincident with bow shock-like features. Archetypes 
of these WR nebulae are those around WR8 and WR18, 
WR35 and WR40 (see Fig. [2] and g]). 

- Al - Mixed WR nebulae. 

These nebulae do not present a clear correspondence be¬ 
tween the optical and IR images. Indeed, they do not show 
a clearly defined IR morphology, with a noticeable lack of 
emission in the WISE W4 band. Examples of these nebulae 
are those around WR 35b and WR 52 (see Fig. O. 

Our classification of the WR nebulae is listed in the last col¬ 
umn of Table [T] We note that sometimes the assignation of a 
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Fig. 1. Examples of WR nebulae with a B type morphology: 
WR 6 (top panels) and WR 16 (bottom panels). Left panels show 
the optical morphology while right panels show the WISE W2 
(blue), W3 (green), and W4 (red) color-composite picture. The 
central WR stars are marked with red lines in the optical images. 
The colour-composite optical picture of WR 6 has been done us¬ 
ing Her (red) and [O iii] (blue) narrow-band images. North is up, 
east to the left. 



Fig. 3. Same as Fig.[T]but for the cases of WR nebulae with a M 
morphology: WR 35b (top panels) and WR 52 (bottom panels). 



Fig. 4. Same as Fig. [T]but for the cases of WR nebulae with a C- 
Disrupted morphology: WR 18 (top panels) and WR 35 (bottom 
panels). 


Fig. 2. Same as Fig. [T] but for the cases of WR nebulae with a 
C-Clumpy morphology: WR8 (top panels) and WR40 (bottom 
panels). The color-composite optical image of WR40 has been 
done using Her (red) and [O iii] (blue) narrow band images. 


morphological type is difficult given the disparity between the 
optical and IR images or due to the ambiguous association be¬ 
tween the star and the nebula. For instance, the nebula around 
WR 116 (Fig. lA.lTl) has been assigned a morphological type C 
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in basis of its optical image, but we reckon that a morphological 
type M cannot be completely excluded by its IR counterpart. 

4.1. Notes on individual objects 

4.1.1. WR35 and its obscured shell 

We report the previously unknown partial shell around WR35. 
This shell, clearly detected in the WISE W4 band (Fig.lU-bottom 
panels), has a size ~2.'0xR5 and opens towards the southeast 
of WR 35. In contrast, the Ha image from the Super COSMOS 
Sky Survey only detects emission from an ionized cloud south¬ 
west of WR 35 which does not coincide spatially with the IR 
shell. A sin iilar case of detection of an obscured shell has been 
claimed by IWachter et al.l (l201ll) for WR 8, although an optica l 
nebulosity was previously reported bv IStock & Barlowl (l20ld) . 
The high frequency of discovery of obscured shells around mas¬ 
sive stars in Spitzer MIPS 24 pm images has led to the sugges- 
tion that such shells may be ubiquitous among these stars (e.g., 
iGvaramadze et al.Ll2010atlWachter et all 120101) . 

4.1.2. A bipolar nebula around WR 85 

The nebula RCW 118 around WR85 has been described as a 
bipolar nebula based on Ha images obtained by iMarston et al.l 
(Il994allbh . This bipolar structure is hinted in the WISE W4 im¬ 
age by the red patches to the northeast and southwest of WR 85 
(Fig. lA.lf] -ng/zr). This IR image is certainly reminiscent of the 
optical ima ge of the nebula NGC (5164-5 around the 0 6 star 
HD 148937 (iBruhweiler et allfl^lDufour et al.Lll988l) . 

4.2. Comparison with previous morphological classifications 

Previous imaging studies of WR nebulae have disclosed a 
rich morphological variety among theni (e.g. , IChu et al.l 119831 
iGruendl ct all |200Q| : I Stock & Barlowl |201Q|) . The most com¬ 
plete classifi cation schem e of WR nebulae is that originall y in¬ 
troduc ed by IChul (Il98ll) and more recently described by Ichul 
(l2003l) . This classification scheme relies on narrow-band opti¬ 
cal images, kinematics, and information on the chemical abun¬ 
dances to define the presence of bubbles and their expansion rate 
and kinematical age. Meanwhile, our classification scheme of 
WR nebulae is supported by the comparison between narrow- 
band optical (mostly Ha) and mid-IR WISE images, the latter 
being very useful to unveil obscured shells and to disentangle the 
nebular emission from background diffuse and ISM emission as 
discussed in Section [3] 

The comparison between these two morphological classifica¬ 
tions of WR nebulae in Table |2] reveals a good agreement. This 
is excellent for the wind-blown bubbles (our !B nebulae and the 
W class), as these are the most easy to spot either in direct im¬ 
ages or kinematical data. Indeed, all our !B nebulae are classified 
as W, and only 4 W objects are classified as C nebulae because 
their mid-IR morphology. 

This consistency is confirmed by the comparison between 
optical im ages in the Ha and [Q ni] emission lines that has been 
shown by iGruendl et al.l (l200Ql) to be very helpful to reveal the 
expansion of shock fronts outside the main nebular shell. For 
example, in the case of RCW 104 around WR75 (Figure lATOb , 
the Ha line image does not show any morphological character¬ 
istic of a round shell, thus implying a classification as a C-type 
nebula. Instead, the [O iii] morphology implies a ^-type. On the 
contrary, the nebu la M 1-67 around WR 124 sho ws no evidence 
of an [O III] front (iFernandez-Martm et al.Ll20l'3l and references 


Table 2. Comparison of WR Nebula Morphology Classification 


WR 

This paper 

Previous works 

References 

6 


WIE 

1,2,3,4 

7 

fB 

W 

1,3,4,5 

8 

C 

E 

4 

16 


W 

4,7 

18 

C 

W 

1,3,4,8 

22 

M 

W 

3,6 

23 


W 

1,3,7 

38 

M 

Ring? 

7 

40 

C 

E 

1,3,4,8 

52 

M 

Rs 

1,3,9 

54 

M 

Ring? 

6 

55 

C 


1,3,9 

68 

M 

Ring 

6 

75 


WIE 

1,3,4,8 

85 

C 

Rs 

1,3,9 

86 

M 


6 

94 

M 


7 

95 

C 


6 

101 

C 

Ring? 

7 

102 


WIE 

4,7 

124 

C 

E 

1,3,4,10 

128 


W/Rs 

3,11 

131 

C 

Hii 

1,3,10 

134 

C 

W 

1,3,12 

136 


WIE 

1,3,4,12 


Refs.- (1) IChul (fl^. (2) IChu et al l (IT^. (3) IChu et all 
(fl^. (4) IStock & Barlo\^ <2010l). (5) llreffers & Chul <1982ah 
(6) |M arston et a 1.1 (Il994al)^(7) IMarston et all (Il994hh, (8) IChul 
(119821). (9) IChu & Treffer j (Il981bh. HOI IChu & Treff^ (Il981ah . 
( 11) IGruendl et all (120001 ). (12) iTreffers & Chul (Il982bl) . 


therein), but its kinematics reve als an expanding shell consist- 
ing of numerous condensations (ISolf & Carsentvtll982h . While 
these findings may raise doubts for the classification of WR neb¬ 
ulae for which [O iii] images and kinematics are not available, 
we note that in all these nebulae the outer expanding shell is 
noticeably much fainter than the shells detected in our ^-type 
nebulae. This points to an intrinsic difference between the two 
morphological types. 

As fo r the M ne bulae, only two have previous definite mor- 
phology (IChuLll99l1) . The nebula around WR 22 is described by 
\C^ (I199¥ as a possible wind-blown bubble, although no evi¬ 
dence for such a bubble is revealed in WISE images (Figure lA^ . 
Similarly, the nebula around WR 52 is categorized as an neb¬ 
ula in basis of kinematical data, but the nebula is rather incon¬ 
spicuous in the IR (Figure [3]-bottom panels). 

4.3. Interpreting the morphology of WR nebulae 

The main morphological characteristics of WR nebulae re¬ 
sult from the interactions of the WR wind and the in¬ 
tense UV radiation with the slow wind ejected previ¬ 
ously by the progenitor. The most basic evolutionary path 
of a WR nebula, as reproduc e d by multiple simulations 


(e.g., 

Garcfa-Segura & Mac LowL 119951: iGarcia-Segura et al.. 

1996c 

JblilFrever et al.. 120031:Ivan Marie et al.l.l2005l: Frever et al.. 

2006: 

lArthuil 2007 allb: Ivan Marie et al.L 120071: Dwarkadas. 

2007; 

iToala & Arthud. l201ll). implies a massive star electing a 


significant fraction of its mass during a RSG or LBV phase 
through a spherical, slow (10-100 km s“^) and dense wind. In the 
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Fig. 5. Spitzer IRS SH and LH spectra of the nebulae around WR7, WR 8, and WR 31a. The spectral responses of the W3 (green) 
and W4 (red) WISE bands normalized to arbitrary values are overplotted on the spectra. The ne bular spectrum of WR 7 corresponds 
to the patch of bright red emission (the WISE W4 band) towards the northwest shown in Figure lATl This spectrum has been plotted 
at two different intensity levels to highlight the low level emission. The contribution of the stellar emission of WR 8 and WR 31a to 
the nebular spectra can be roughly reproduced by a black-body model (dashed line). 


hnal stage of stellar evolution, during the WR phase, this CSM 
will be swept up by the WR fast stellar wind (Voo >10^ km s“^), 
creating at early times a WR bubble (^-type nebulae). Later on, 
the WR bubble will be disrupted due to Rayleigh-Taylor and thin 
shell instabilities (C-type nebulae) and the resulting hlaments 
and clumps will later mix with the ISM (Al-type nebulae). 

Additional processes may alter this general scheme. For in¬ 
stance, the proper motion of the star or a density gradient in 
the ISM may compress and brighten the nebula as it interacts 
with its surroundings, resulting in incomplete shells and arcs. 
Stellar rota tion can also introduce changes in the s tellar evo¬ 
lution fe.g.. iMevnet & Maedeii 120031: lEkstrom et al.L|20T^ and 
produce, for instance, anisotropic stellar winds. In the case of 
non-spherical (e.g., bipolar or clumpy) ejections of material seen 
in some LBV nebulae (e.g., M N13 and MN79; iHumphre"^ 
|20 1 (A [^aramadze et al .L 12010^ . the wind-wind interaction will 
give place to the immediate formation of clumps or hlaments. 


forming directly a C-type nebulae with no B-type nebula in be¬ 
tween. The large-scale structure of the ISM may also alter the 
morphology of WR nebulae. 

Since WR nebulae can form through different morphological 
sequences, we should not expect a tight correlation between neb¬ 
ular morphology and t he star evolutionary stage ini plied by its 
WR spectral type fe. g.. iMof^ 1 1 9951: iMaedeii 1 1 9971) . Certainly, 
Al-type nebulae can be expected preferentially around evolved 
WR stars, whereas the central stars of ^ and C-type nebulae can 
be ex pected to have late spectral types (IChuL 1 19811: IChu et al.l 
Il983h . This correspondence is further complicated by the time- 
lapse spent by a WR star on the different WR sub-phases which 
depends strongly on different factors (initial mass, rotati o n, ini¬ 
tial metahicity; e.g ., iMaederl Il99ll: iMevnet & Maeded. 120051: 
iGeorgv et al.L l20l3) . A star may even miss a WR sub-phase; 
for example, the stellar e volution models at solar metahicity of 
iMevnet & Maeded (l2003l) predict that a rotating 25 Mq star will 
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Toala et aL: WISE morphological study of WR nebulae 



WNL WNE WC 


Fig. 6. Distribution of the different morphological types of WR 
nebulae depending on the spectral types of their WR central 
stars. 


not enter the WNE and WC stages. Therefore, for those particu¬ 
lar models, any type of WR nebulae can be found around WNL 
stars because the long duration of this phase or even the lack of 
the WNE stage. On the other hand, for the case of non-rotating 
models, which predict comparable durations of the WNL and 
WNE phases, we may expect ^-type nebulae around these two 
type of stars. 

To better quantify these relationships, we show in Eigure[6| 
the distributions of nebular morphology and stellar spectral type 
among the WR nebulae in our sample. Live out of the nine B 
nebulae have WNE stars and two others have WNL spectral 
types, whereas only two are early WC stars. The WNL stars 
WR 16 (WN8) and WR31a (WNll) are surrounded by small¬ 
sized bubbles (2.8 and 1.2 pc in radius, respectively), and the 
WC WR23 (WC6) and WR 102 (W02) stars by large bubbles 
(16 and 4 pc in radius, respectively). Nebulae with C-type mor¬ 
phology harbor mostly late WN7-WN8 stars (6 out of 12) or 
WNE stars (4 out of 12), with only two WC stars (WR95 and 
WR 101). On the contrary, M nebulae have a significant fraction 
of early WC stars (6 out of 10) and a very small number of WNL 
stars (1 out of 10). This confirms previous suggestions that WR 
bubbles are associated to late WN stars and C-type to WC stars. 
Most C nebulae had been previously classified either as W or 
nebulae, but we notice some contamination of and neb¬ 
ulae, and even objects of uncertain morphology. This indicates 
that C nebulae have a variety of kinematical and chemical prop¬ 
erties. This is also the case for their WR central stars which can 
be linked to different evolutionary paths leading to a C-nebula 
morphology. 

5. Comparison with previous mid-IR studies 

The extinction by dust in the Galactic Plane makes difficult the 
search for WR nebulae at optical wavelengths. The detection 
of the mid-IR emission (thermal dust, molecules, PAHs, ...) of 
shells around evolved stars may provide an alternative method 
for the direct identification of these nebulae. Mid-IR observa- 
tions have indeed allowed the detection of yo ung PNe (e.g., 
iMorris et"^ 120061: lEesen & MilisavljevicLl201Qh and LBV neb¬ 


ulae fe.g.. iGvaramadze et al.L I20l65 IWachter et al.L l201Ql) . Not 
unexpectedly, Spitzer MIPS 24 pm data have been used to dis¬ 
cover new WR stars by the direct detection of mid-IR emission 
from it s circumstellar nebula (e.g., WR 121b. lGvaramadze et al.L 
l2010ch . Two outstanding examples are the identification of the 
interacting neb ulae around the WN9h typ e stars, WR 120bb and 
WR 120bc, by iBurgemeister et ^ (l2Q13l) and the identification 
of a WR nebula i n the Large Magellanic Cloud reported by 
iGvaramadze et al.l (l2014l) . Of course, the direct identification of 
circumstellar nebulae around WR stars in the mid-IR 24 pm 
band requires a mandatory spectroscopic characterization of the 
nebular emission, but it is important to emphasize the simplicity 
and potential of the mid-IR method. 

In this paper, we have taken advantage of the similarities 
between the spectral responses of the Spitzer MIPS 24 pm 
and WISE W4 22 pm bands to investigate the circumstellar 
medium around a significative sample of known Galactic WR 
stars. By identifying the mid-IR counterparts of these WR neb¬ 
ulae, we have assessed their true extent, disentangling the op¬ 
tical emission of the circumstellar shells from that of the ISM 
along the line of sight. New mid-IR shells (e.g., WR35, Lig.lH- 
bottom panels) have also been discovered. It is worth mentioning 
here the notable case of NGC 3199, the n ebula around WR 18 
(see Pig.|4]-top panels). [Stock et al] (|2Q1 ll) derived the chemical 
abundances for a bright Ha nebular clump and concluded they 
were consistent with those of the Galactic H ii region M 17 on 
which NGC 3199 is projected. Accordingly, they proposed that 
NGC 3199 consists mainly of swept-up ISM material. A careful 
registration of the location of the ~15" VLT/UVES slit used by 
IStock et akl (l2Qlll) on the WISE images shows that this small slit 
does not include the emission detected in the W4 band (i.e., neb¬ 
ular material), but it is located on a spot of bright emission in 
the W3 band (i.e., ISM material). These findings cast doubt on 
the ISM abundances implied by those authors fo r NGC 3199 and 
rather support those reported bv lMarstc^ (1200 ll) . 

6. Summary and conclusions 

We have examined WISE IR images of a sample of 31 WR stars 
to study their nebular morphologies in mid-IR bands. We have 
then compared the emission in the WISE images with optical Ha 
and (in some cases) [O iii] emission line images. The variety of 
morphologies of WR nebulae can be classified as bubble ^-type 
nebulae, clumpy/disrupted C-type nebulae, and mixed Al-type 
nebulae. 

We have used Spitzer IRS spectra of some WR nebulae to 
investigate the nature of the IR emission of WR nebulae in the 
different WISE bands (Toala et al. in preparation). According to 
previous IR studies of evolved massive stars, we find that the 
emission in the WISE W4 band of tB and C-type nebulae most 
likely traces dust associated to the WR nebula. On the other 
hand, the WISE W3 band is mainly tracing the emission of ma¬ 
terial along the line of sight of the nebula, which results in di¬ 
minished emission at optical wavelengths. We emphasize that 
the acquisition of mid-IR images of WR nebulae, especially for 
wavelengths above 20 pm, is very useful to disentangle the emis¬ 
sion from the WR nebula from that of the ISM. This provides 
a less demanding observational approach for the identification 
and study of WR nebulae than those based on optical images 
and spectra. The discovery of an obscured shell around WR 35 
which is detected only at 22 pm in the W4 band WISE images or 
the capability to distinguish the nebular material of NGC 3199 
around WR 18 from that of the Galactic H ii region M 17 clearly 
illustrate the advantages of mid-IR observations of WR nebulae. 
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Toala et aL: WISE morphological study of WR nebulae 


We find a loose correlation between nebular morphology and 
stellar spectral type, also claimed in the past, but we note that 
this correlation is complicated by the different evolutionary se¬ 
quences of WR nebulae and the dependence of the stellar evo¬ 
lution in these phases with initial mass and composition, and 
stellar rotation. 
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Appendix A: Additional Figures 

In this appendix we collect the figures of the rest of the WR 
nebulae studied in this paper not presented in the main text. The 
nebulae are those around the WR stars WR7, WR22, WR23, 
WR30, WR31a, WR38, WR54, WR55, WR68, WR75, WR85, 
WR86, WR94, WR95, WR 101, WR 102, WR 113, WR 116, WR 124, 
WR 128, WR 131, WR 134, and WR 136. 



Toala et aL: WISE morphological study of WR nebulae 



Fig. A.l. Optical (left) and mid-IR WISE (right) images of WR 7 
(NGC2359). See Table 1 for details of the optical image. The 
central WR star is marked with red lines in the left panel. North 
is up, East to the left. 



Fig. A.5. Same as Fig. lA.ll for WR 31a. 



Fig. A.6. Same as Fig. lA.ll for WR 38. 



Fig. A.3. Same as Fig. [Allfor WR 23. 





Fig. A.8. Same as Fig. [Allfor WR 55 (RCW 78). 



Fig. A.4. Same as Fig. lA.ll for WR 30 (Anon). 
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Fig, A.IO. Same as Fig.Ofor WR75 (RCW 104). 


Fig. A.14. Same as Fig. lA.ll for WR95. 



Fig, A,ll, Same as Fig.Ofor WR85 (RCW 118). Fig, A,15. Same as Fig.Ofor WR 101 (Anon). 



Fig, A,12, Same as Fig.Ofor WR 86 (RCW 130). 


Fig, A,16. Same as Fig.Ofor WR 102 (G2.4+1.4). 
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Fig. A.18. Same as Fig. lA.ll for WR 116 (Anon). 


Fig. A.22. Same as Fig. lA.ll for WR 134 (Anon). 




Fig. A.19. Same as Fig.Ofor WR 124 (Ml-67) 


Fig. A.23. Same as Fig.Ofor WR 136 (NGC 6888). 



Fig. A.20. Same as Fig. lA.ll for WR 128 (Anon). 
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